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MECHANICS  OF  COMPOSITE  MATERIALS^ 

by 

G.  A.  Hegemier^ 
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ABSTRACT 

This  report  summarizes  the  results  of  a theoretical  study 
of  wave  propagation,  diffusion,  and  damage  accumulation  in 
composite  materials. 


Sponsored  by  the  United  States  Air  Force  Office  of  Scientific  Research, 
at  the  Univ'ersity  of  California,  San  Diego. 


Professor,  Department  of  Applied  Mechanics  and  Engineering 
Sciences,  University  of  California,  San  Diego,  La  Jolla,  California. 


MECHANICS  OF  COMPOSITE  MATERIALS 


1.  RESEARCH  OBJECTIVES 

Research  under  the  subject  Contract  represents  a continuing  effort^ 
to  construct  viable  nonlinear  micro-mechanical  continuum  models  to  pre- 
dict the  thermomechanical  response  of  advanced  composite  materials.  The 
thermomechanical  processes  under  study  include: 

. Linear  and  nonlinear  wave  propagation 

. Linear  and  nonlinear  diffusion  (thermal  and  hygrothe rmal) 

. Mechanisms  and  accumulation  of  damage. 

The  thermomechanical  measures  include: 

. Global  stress,  deformation,  temperature,  internal  energy  fields 
. Local  stress,  deformation,  temperature,  internal  energy  fields 
. Re  s idua I prope  r ties  of  damaged  composites. 

2.  RESEARCH  PROGRESS 

2.  1 Discus  s ion 

Excellent  progress  has  been  made  to  date  in  modeling  both  laminated 
and  fibrous  composites  as  homogeneous  continua  with  microstructure.  The 
resulting  theory  is  generally  referred  to  as  a "Theory  of  Interacting  Continua" 
or  a "Mixture  Theory  with  Microstructure".  In  contrast  to  classical  theories 

■^Previous  work  was  conducted  under  AF-AFOSR  Grant  70-1975  and  AF-AFOSR 
Contract  75-2870. 
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with  similar  titles,  the  governing  equations  are  completely  specified  by  a 


knowledge  of  the  initial  microstructural  geometry  of  the  composite,  the 
component  constitutive  relations,  and  the  component  interface  physics. 

In  further  contrast  to  classical  theories,  the  current  model  provides  infor- 
mation on  stress,  displacement,  strain,  temperature,  and  internal  energy 
fields  within  the  microcomponents  as  well  as  global  measures  of  these 
quantities.  The  present  theory  also  accounts  for  debonding  at  the  micro- 
component interfaces  - a basic  damage  mode.  Models  incorporating  crack 
distributions  throughout  the  constituents  are  under  study. 

The  most  general  mathematical  model  considered  to  date  is  a non- 
linear, anisotropic,  the rmomechanical  mixture  theory  with  microstructure. 
Particular,  simplified  forms  of  this  model  have  been  developed  to  cover 
special  cases  of  composite  geometry,  material  constitutive  relations,  and 
physical  processes  or  problems.  In  all  cases,  however,  theoretical 
construction  is  based  upon  one  of  two  basic  methods  which  are  in  turn 
based  upon  an  asymptotic  scheme  in  which  dominant  signal  wavelengths 
associated  with  the  physical  process  are  assumed  to  be  large  in  comparison 
with  typical  composite  microdimensions.  An  explanation  of  the  foregoing 
construction  methods,  as  well  as  a number  of  interesting  examples,  can 
be  found  in  Appendix  II. 

The  composites  studied  using  the  above  asymptotic  technique  range 
from  elementary  laminated  to  complex  three-dimensional  advanced 
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composites.  In  addition  to  certain  model  composites,  materials  have 


included  quartz  phenolic,  carbon  phenolic,  carbon  carbon,  and  graphite 
epoxy.  Representativ'e  problem  areas  include  nose  tips,  heat  shields, 
engine  blades,  nozzles,  and  hygrothermal  environments. 

Cases  completed  and  in  publication  thus  far  include: 

. A general  theory  for  wave  propagation  in  linear  elastic  and 
linear- vis  coelastic  laminated  composites  with  isotropic  laminae; 

. A general  theory  for  linear  thermal  diffusion  in  laminated 
composites  with  isotropic  laminae; 

. A "first  order"  theory  for  linear  wave  propagation  in  anisotropic 
laminates  with  oblique  layup; 

. "First  order"  theories  for  wave  propagation,  diffusion,  and 
debonding  in  unidirectional  fibrous  composites  with  cylindrical  fibers  of 
arbitrary  cross  section; 

. "First  order'  nonlinear,  anisotropic  and  isotropic,  thermo- 

I 

dynamic  theories  for  wave  propagation,  diffusion,  and  damage  acciunulation  j 

due  to  debonding  in  both  laminated  and  fibrous  composites  of  specific  types. 

In  the  general  linear  cases  above,  a hierarchy  of  models  is  system-  | 

atically  defined  by  the  order  of  the  truncation  of  the  asymptotic  sequence  ! 

involved  in  the  theoretical  construction  process.  Retention  of  all  terms 
in  this  sequence  leads  to  a formally  exact  theory.  Retention  of  only  the 
zeroth-ordor  plus  first  order  terms  results  in  a so-called  "first  order" 
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theory  which  can  be  cast  in  the  form  of  a nrixture  theory.  The  latter 
usually  constitutes  an  adequate  micromechanical  model  of  a composite 
for  engineering  purposes.  All  nonlinear  models  are  "first  order"  theories. 

The  first  order  models  for  laminates  and  fibrous  composites 
incorporate  the  effects  of  arbitrary  fiber  geometry,  anisotropy,  and  inter- 
face debonding.  In  addition  to  diffusion,  a thermodynamic  theory  of  elastic- 
plastic  wave  propagation  under  finite  deformation  has  been  developed.  The 
constituti ves  laws  of  the  microcomponents  currently  in  publication  relate 
the  mean  stress  (trace  of  the  stress  tensor)  to  density  and  internal  energy 
by  a Mie-Griineisen  caloric  equation  of  state  and  the  stress  deviator  ten- 
sor to  the  rate-of-deformation  tensor  by  an  elastic-perfectly  plastic 
relation  of  the  von  Mises  type;  extension  to  general  anisotropic  plasticity 
utilizing  a plastic  potential  and  an  associated  flow  ride  is  under  develop- 
ment. Constituent  debonding  and  component- inte r face  physics  is  governed 
by  a Coulomb  frictional- type  relation;  the  latter  constitutes  the  present 
m.easure  of  damage";  extension  of  this  measure  to  include  microcracking 
throughout  the  components  is,  as  was  noted  previously,  under  study. 

In  the  case  of  linear  elastic  wave  propagation,  model  accuracy 
has  been  demonstrated  by  comparison  of  phase  velocity  spectra  with  exact 
and/or  experimental  data.  Accuracy  superior  to  existing  theories  has 
been  observed,  as  well  as  first  order  model  applicability  down  to  wave- 
lengths on  the  order  of  typical  composite  microdimensions  (when  the  energy 
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is  partioned  priniarily  in  the  first  mode).  In  addition,  a number  of  transient 
pulse  cases  have  been  treated;  these  exhibit  good  correlation  with  exact  and/or 
experimental  data.  A sampling  of  typical  transient  pulse  results  for  both 
laminated  and  unidirectional  fibrous  composites  is  illustrated  in  Figs.  1-4 
of  Appendix  1.  Typical  transient  wave  propagation  results  for  a complex, 
advanced,  three-dimensional  composite  are  given  in  Figs.  8-10  of  Appendix  I. 
Complete  details,  including  extensive  phase- velocity- spe ctra- cal culations , 
can  be  found  in  the  appropriate  references  listed  under  Section  2.  3. 

In  the  case  of  large  deformations  and  elastic-plastic  wave  propaga- 
tion, theoretical  transient  pulse  - results  for  propagation  parallel  to  the 
fibers  of  a unidirectional  fibrous  composite  have  been  compared  with  data 
from  a 2 D- I.agrangian  finite  difference-code  (CRAM);  the  latter  is  essentially 
exact.  As  can  be  observed  from  Appendix  I,  Figs.  11-14,  the  agreement 
between  "exact"  and  approximate  calculations  was  found  to  be  excellent. 

Model  accuracy  in  the  case  of  interface  debonding  has,  to-date, 
been  assessed  by  comparison  with  experimental  data  on  delaminated 
plates  composed  of  alternating  layers  of  Polymethyl  Methacrylate  (PMMA, 
Rohm  and  I-Iaas  Type  A),  and  6061 -To  aluminum.  The  laminae  of  the 
composite  were  oriented  perpendicular  to  the  impact  plane,  and  struck 
by  a projectile  fired  from  a 10  cm  bore  light  gas  gun.  On  the  basis  of 
averaged  particle  velocity  measured  at  the  rear  face  of  the  specimen, 
the  model  hypothesized  yields  excellent  correlation  with  the  experimental 
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results.  The  test  setup  and  typical  results  are  illustrated  in  Figs.  15-18 
of  Appendix  I, 


In  addition  to  the  above,  studies  concerning  the  accuracy  of  the 
thermal-portion  of  the  continuum  models  were  undertaken.  In  particular, 
diffusion-type  processes  in  a laminated  composite  with  periodic  micro- 
structure wore  examined.  Solutions  for  the  lowest-order  niodels  were 
compared  with  'exact'  results  from  a finite  difference  code.  P'or  most 
cases  the  lowest-order  "eftective  conductivity"  theory  was  found  to  yield 
good  results.  For  exceptional  problems  requiring  a higher-order  theory, 
a modified  version  of  the  lowest-order  theory  was  found  to  yield  excellent 
correlation  between  exact  and  approximate  solutions.  For  many  problems 
of  the  diffus  ion-type,  these  elementary  equations  may  offer  an  attractive 
alternative  to  other  means  for  obtaining  solutions.  Typical  results  for 
the  case  of  heat  propagating  norm.al  to  the  layers  of  the  laminated  com- 
posite of  Fig.  1,  Appendix  1,  are  illustrated  in  Figs.  19  and  20  of 
Appendix  I.  In  each  case  a uniform  temperature  T is  applied  to  y = 0 
from  t = 0 and  t - t^  while  the  boundary  temperature  is  zero  for  all  other 

time.  (Here  A =h  + h , p c = heat  capacity  of  or-material,  k - 

1 2 o O a 

thermal  diffusi vity).  It  is  emphasized  here  that,  although  F'igs.  19,  20 
refer  to  a one-dimons  ional  example,  the  theory  is  valid  for  three- 
dimensional  problems.  The  one-dimensional  case  was  selected  mainly 
because  the  most  severe  test  of  the  theory  occurs  when  heat  flow  is 
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normal  to  the  laminates  (since  material  properties  are  discontinuous  in 
this  direction). 

finally,  it  is  noted  that  a major  development  in  theoretical  modeling 
was  effected  during  the  previous  research  period.  This  development  allows 
treatment  of  arbitrary  fiber  geometry  and  arrays.  Basically  the  task  of 
constructing  mixture  theories  has  been  partioned  into  two  isolated  steps. 

The  first  step  involves  the  determination  of  the  gene ral  mixture  equations 
for  a particular  problem  class;  these  equations  involve  coefficients  which 
must  be  determined  from  the  detailed  microstructural  geometry,  inter- 
face data,  and  con3tituti\e  relations.  The  second  step  defines  the  con- 
stants via  a quasi-static  ' rn  ic  robounda  ry  value  problem"  (MV’P). 

Resolution  of  the  is  carried  out  by  a variational-based  finite  element 

n'.et'nod.  The  procedure  has  been  applied  to  diffusion  and  wave  propaga- 
tion in  unidirectional  composites.  Here  important  information  has  been 
obtained  regarding  the  accuracy  of  the  concentric  cylinders  approximation 
used  frequently  in  practice  for  circular  fibers  in  a hexagonal  array  and 
for  rectangular  fibers  in  a similar  array.  Extension  of  the  procedure  to 
arbitrary  three  dimensional  geometries  is  under  development. 

A discussion  of  the  above  subject  matter  can  be  found  in  the  survey 
paper  listed  under  Section  2.  3.  F'or  the  readers  convenience,  the  latter 
IS  included  as  Appendix  II  of  this  report. 
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2.  2 Research  Applications 

The  fundamental  concepts  embodied  in  the  aforementioned  wave 
propagation  studies  have  been  successfully  employed  ir?  an  applied  research 
effort  to  model  1-D  (FMI  COMRAD)  and  AVCO  3DQP  quartz  pheonlic 
composites.  These  niodels,  and  the  resulting  numerical  code  (TINC)  are 
presently  in  practical  use  by  the  engineering  profession.  Typical  results 
for  3 DQP  are  illustrated  in  Figs.  5 - 10  of  Appendix  I. 
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APPENDIX  I 


Compend Liuii  of  Typical  Theoretical  Results 


Preseated  iti  this  Appendix  is  a san\pling  of  typical  transient  pulse 
and  diffusion  calculations  carried  out  using  the  Theory  of  Interacting 
Continua  (TINC).  The  cases  discussed  include  laminated  and  unidirectional 
fibrous  composites  as  \^cll  as  a complex  3 -dimens ional  layup.  Results  of 
the  TINC  theory  for  each  example  are  compared  to  experimental  data. 


l.aminated  Composites 

The  geometry  of  the  laminated  composites  studies  is  illustrated 
in  I'ig.  1.  One  problem  of  interest  is  that  of  symnaetric  (P)  wave 
propagation  parallel  to  the  layers,  with  material  behavior  restricted  to 
the  linear  clastic  range.  A comparison  of  TINC  and  experimental  results 
is  given  in  P'ig.  2.  The  experimental  results  were  reported  by  Whittier 
and  Peck  [1  j for  specimens  composed  of  Thornel  (high  modulus  carbon) 
fibers  reinforcing  a carbon  phenolic  matrix.  Specimens  of  1/4-inch 
thickness  were  subjected  to  a uniform  pressure  at  the  left  boundary,  with 
a step  function  time- dependence  induced  by  a gas  dynamic  shock  wave  of 
about  70  psi.  The  numerical  calculations  were  initiated  by  impacting  a 
step  function  velocity  of  7.786  cm/sec  to  both  constituents  at  the  boundary. 
While  this  condition  does  not  correspond  precisely  to  the  experiment,  it 
was  felt  that  the  error  introduci^d  would  be  negligible  when  far  removed 
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from  the  front  surface.  Figure  2 depicts  the  comparison  between 
experiment  and  theoretical  results  for  averaged  rear  surface  velocity, 
normalized  on  the  boundary  input  velocity.  Since  absolute  times  were 
not  measured  in  the  tests,  theoretical  and  experimental  results  were 
matched  at  their  respective  first  peak  arrivals.  In  the  figure,  results 
from  the  full  TINC  equations  solved  numerically  are  denoted  "Continuum 
mixture  theory."  The  results  denoted  "Simplified  theory"  are  for  an 
elementary  form,  of  the  TINC  theory  admitting  a closed-form,  solution. 

Unidirectional  Fibrous  Composite 

The  TINC  theory  has  been  applied  to  the  problem  of  wave  prop- 
agation parallel  to  the  fibers  of  a unidirectionally  reinforced  composite 
with  hexagonal  array.  The  hexagonal  array  is  appro.ximated  by  concentric, 
li.nearly  elastic  cylinders  as  illustrated  in  Fig.  3. 

In  Fig.  4,  transient  pulse  data  predicted  by  TINC  is  compared 
with  experimental  data  on  a unidirectional  quartz  phenolic  fibrous 
composite.  The  experiments  were  subjected  to  a 70  psi  step  function 
in  pressure  via  a shock  tube.  The  figure  depicts  the  comparison  between 
the  experimental  and  TINC  theory  code  predictions  of  a r ea -a ve raged  rear 
surface  velocity.  Results  of  the  simplified  TINC  theory  arc  also  por- 
trayed. Once  again  the  experimental  and  theoretical  results  are  matched 
at  their  first  peak  arrivals  since  absolute  arrival  times  were  not  measured 
experimentally.  As  was  the  case  in  Fig.  2,  the  times  shown  arc  those 


predicted  by  the  theory. 


Three-dimensional  Layups 

Predictions  of  transient  pulse  behavior  from  the  TINC  theory 
have  been  compared  to  experimental  results  for  a 3-dimensional  quarts 
phenolic  composite  (3DQP)  manufactured  by  AVCO.  A representative 
selection  of  the  results  are  given  here. 

A "block"  3DQP  specimen  (Fig.  5)  was  subjected  to  shock  tube 
tests  by  the  Aerospace  Corporation.  The  test  procedure  is  the  same 
as  that  used  for  laminated  and  fibrous  composites  and  is  described  above. 
Once  again  the  theoretical  and  experimental  boundary  conditions  are  not 
identical,  but  this  decrepancy  is  expected  to  introduce  little  error.  The 
averaged  rear  surface  velocity  measured  by  experiment  is  compared  to 
that  predicted  by  the  TIXC  theory  in  Fig.  6.  Also  shown  are  results 
predicted  by  the  elementary  "Flead  of  the  Pulse"  approximate  [2]. 

The  experimental  data  has  been  forced  to  coincide  with  the  theoretical  at 
^ ^ ^ Q ~ ^ ^ ^ since  absolute  arrival  times  we  re  not  measured. 

Another  set  of  relatively  low  pressure  wave  propagation  experiments 
w-ere  performed  at  the  Air  Force  Weapons  Laboratory  using  a light  gas 
gun  and  "midspace"  3DQP  specimens,  shown  in  Fig.  7.  In  this  test,  three 
specimens  of  different  thicknesses  were  mounted  on  a single  target  and 
impacted  with  a ,0  1 -inch  mylar  flyer.  The  input  stress  profiles  for  the 
three  specimens  are  illustrated  in  Fig.  8.  For  the  TINC  calculations. 
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input  velocity  time-histories  corresponding  to  the  stress-time  data  of 
Fig.  8 were  used.  A comparison  of  experimental  and  theoretical  results 
for  the  thickest  of  these  specimens  is  shown  in  Fig.  9.  Since  absolute 
arrival  times  were  not  measured,  the  two  curves  are  forced  to  coincide 
at  their  peak  pressures.  The  agreement  here  is  excellent,  especially 
in  view  of  the  fact  that  the  wave  is  measured  only  2-1/2  microdimensions 
from  the  front  surface.  Also,  with  the  lOkb  input  stress,  this  is  a non- 
r?  1 rr  ti  on, 

The  AFWL  also  performed  a series  of  high-pressure  gas  gun 
experiments  on  ''midspace''  3DQP  specimens  (Fig.  7).  Experimental 
and  theoretical  data  are  compared  in  Fig.  10.  The  observed  agreement 
indicates  that  the  thermodynamic  and  nonlinear  constitutive  behavior 
incorporated  into  the  TINC  theoretical  model  are  adequate.  Note  that 
here  the  rear  surface  velocity  is  measured  only  1-1/2  microdimensions 
from  the  front  surface.  Thus,  as  before,  the  TINC  theory,  though 
derived  as  an  "outer"  solution,  yields  surprisingly  good  results  even 
at  locations  quite  near  the  boundary. 

Comparison  of  TINC  and  2-D  Finite  Difference  Code 

It  is  of  interest  to  compare  the  results  from  a TINC  calculation 
to  those  of  a typical  two-dimensional  Lagrangian  finite  difference  code 
which  is  essentially  exact.  Such  a code  (CRAM)  is  available  at  Systems, 
Science  and  Software,  Inc.  The  two  codes  were  compared  for  the 
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concentric  cylinder  geometry  of  Fig.  11.  Computer  costs  for  CRAM  were 
approximately  500  times  those  for  running  the  TINC  code.  Several  typical 
results  are  illustrated  in  Figs.  12-14  for  a square  wave  input  velocity 
(shown  as  problem  jfZ  on  Fig.  11).  Once  again,  the  calculations  were 
carried  out  within  a few  microdimensions  of  the  boundary. 

Debonding  of  Constituent  Interfaces 

Model  accuracy  in  the  case  of  interface  debonding  or  cracking 
has,  to  date,  been  assessed  by  comparison  with  experimental  data 
dealing  with  delaminated  plates  composed  of  alternating  layers  of 
Polymethyl  Methacrylate  (PMMA,  Rohm  and  Haas  Type  A),  and  6061-T6 
aluminum.  The  laminae  of  the  composite  were  oriented  perpendicular 
to  the  impact  plane,  and  struck  by  a projectile  fired  from  a 10  cm  bore 
light  gas  gun.  On  the  basis  of  averaged  particle  velocity  measured  at 
the  rear  face  of  the  specimen,  the  model  hypothesized  yields  excellent 
correlation  with  the  experimental  results. 

The  test  setup  employed  is  illustrated  in  Fig.  15.  Typical 
correlation  between  theoretical  and  experimental  results  is  shown  in 
Figs.  16-18.  Complete  details  can  be  found  in  reference  [13^ 

Dif  fus  ion 

The  laminated  geometry  of  Fig.  1 was  selected  to  evaluate  the 
continuum  model  construction  procedure  for  diffusion  problems  (a 
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complete  discussion  of  this  work  can  be  found  in  reference  [4]).  The 


case  of  heat  propagation  normal  (in  the  y-direction)  to  the  layers  was 
selected  for  a comparison  of  exact  and  continuum  theory  solutions  since 
it  is  known  that  this  case  represents  the  most  severe  test  of  the  latter. 
Figure  19  represents  a case  of  similar  material  properties  whereas 
Fig.  20  illustrates  the  comparison  for  widely  differing  material  prop- 
erties. Continuum  theory-results  were  obtained  via  closed-form 
■?n)ininns  for  a snua  re-wave  temperature  input  at  the  boundary  y = 0. 

Exact"  results  were  obtained  using  a finite  difference  code.  The 
calculations  illustrate  the  ability  of  the  first-order  theory  to  model 
microstructure  details  and  closely  match  the  exact  data. 
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Conparison  of  TINC  and  experimental  rear  surface 
velc.city  results  from  quartz  fiber  reinforced  phenolic 


Fi^.  13.  Comparison  of  'I  INC  and  avorago  CRAM  particle 
- "velocity  histories  for  sheath. 
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Fig.  17.  Comparison  of  theory  with  experiment  No.  6. 
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Fit’.  19.  Temperature  distribution  at  r - I,  1.5  and  2 for  similar 
material  properties  and  sniall  duration  of  input  pulse  t . 
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MIXTURE  THEORIES  WITH  MICROSTRUCTURE  FOR  WAVE 
PROPAGATION  AND  DIFFUSION  IN  COMPOSITE  MATERIALS 


Two  procedures  for  constructing  mixture  theories  with  micro- 
structure are  discussed  for  wave  propagation  and  thermal  diffusion 
processes  in  binary  laminated  and  fibrous  composites.  Each  method 
proceeds  from  the  composite  microstructure  and  is  based  upon  an 
asymptotic  scheme  wherein  the  ratio  of  trans verse-to-longitudinal 
characteristic  times  associated  with  a physical  process  is  assumed  to  be 
small.  The  methods  retain  considerable  information  regarding  the 
m.echanical  and/or  diffusive  fields  within  the  microcomponents. 

Theoretical  results  are  compared  with  both  experimental  and  exact  data 
in  an  effort  to  evaluate  model  accuracy.  The  presentation  includes  treat- 
ment of  geometric,  constitutive,  and  interface  nonlinearities. 

l.  INTRODUCnON 

A considerable  number  of  continuum  theories  have  been  proposed 
for  thermomechanical  processes  in  composites.  The  continuum  theory 
of  mixtures  represents  one  of  the  more  Important  and  successful  theoreti- 
cal descriptions  of  such  multiphase  materials.  According  to  this  concept, 
the  composite  constituents  are  modeled,  at  each  instant  of  time,  as  super- 
posed contlnua  in  space.  Each  continuum  is  allowed  to  undergo  individ- 
ual de form.ations.  The  microstructure  of  an  actual  composite  is  then 
simulated  by  specifying  the  interactions  between  the  continua. 

Mixture  theories  for  multiphase  materials  have  been  developed 
from  both  macrostructural  and  m ic ros tructural  viewpoints.  The  vast 
majority  of  activity,  however,  has  been  directed  toward  the  former, 
which  concerns  the  formulation  of  a general  theoretical  framework  for 
classes  of  m.ixtures  based  upon  certain  phenomenological  postulates  that 
avoid  detailed  microstructural  considerations.  Excellent  accounts  of 
this  approach  can  be  found  in  [ 1-3]. 

The  development  of  mixture  theories  from  the  microstructural 
viewpoint  represents  a m.o  re  difficult  task.  The  reward,  however,  is  a 

m. odel  - which  shall  be  called  a mixture  theory  with  microstructure  - 
that  is  corr.pletely  determined  by  the  geometry  and  constitutive  relations 
of  the  individual  components  and  which  provides,  to  a certain  degree  of 
accuracy,  information  concerning  microstructure  fields.  0‘  course, 
the  price  one  must  pay  for  such  detail  is  often  a case-by-case  study 
rather  than  a general  theory.  Nevertheless,  it  should  be  recognized 
that,  while  general  macrostructural  - based  theories  may  appear  more 
elegant,  such  theories  may  not  be  practical  for  many  applications  due  to 
the  unreasonable  burden  placed  upon  the  experimentalist  to  determine  a 
multitude  of  constants  and/or  functions. 
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In  this  presentation,  focus  is  placed  upon  several  recent  works  b/ 
the  author  and  co-workers  concerning  the  construction  of  mixture  theories 
^vith  microstructure  for  laminated  and  fibrous  composites.  The  discus- 
sion is  limited  to  thermomechanical  processes  in  composites  with  initial 
periodic  microstructure.  Particular  emphasis  is  placed  upon  wave 
propagation  and  thermal  diffusion. 

The  presentation  is  divided  into  five  sections.  The  methodology- 
associated  with  two  theoretical  construction  techniques  is  given  in 
Section  2 along  with  a brief  survey  of  applications.  These  techniques 
are  demonstrated  via  two  examples  in  Sections  3 and  4.  Finally,  in  an 
effort  to  illusidate  the  applicability  and  utility  of  the  techniques,  tv/o 
applications  to  nonlinear  problems  are  discussed  in  some  detail  in 
Section  5;  these  include  comparisons  between  mixture  theory,  exact  and/ 
or  experimental  results. 

2.  METHODOLOGY  AND  SURVEY 

A variety  of  construction  procedures  may  be  used  to  generate 
mixture  theories  with  microstructure.  Two  such  procedures  are  selected 
dia\.uooicn  hci'ci.n  in  vic'.v  of  tlioir  observed  success  in  describing 
thermomechanical  processes  in  composite  materials. 

The  first  procedure,  which  is  Illustrated  in  Section  3 by  an 
elementary  example,  is  the  Dif fe r ential- dif fe r ence  Method.  This 
approach  has  been  used  to  develop  continuum  models  of  wave  propagation 
in  linear  elastic  r4-9,2D]  and  linear  viscoelastic  [9-10]  laminated  com- 
posites, wave  propagation  in  linear  elastic  fibrous  composites  [7-11], 
w’ave  propagation  and  debonding  in  laminated  composites  [ 12],  and 
thermal  diffusion  in  laminated  composites  [13,  14].  The  method  consists 
of  eight  basic  construction  steps  which  may  be  summarized  as  follows: 

1)  the  field  quantities  (all  dependent  variables)  of  each  composite 
constituent  are  expanded  in  a suitable  spatial  series  (this  expansion, 
which  amounts  to  discretization  of  the  component  domains,  may  range 
from  a Taylor  series  for  sufficiently  simple  geometries  to  topological 
multiplexes  'IS],  e.  g.  , finite  elements,  for  complex  geometries); 

2)  recurrence  relations  for  the  expansion  coefficient-functions  are  ob- 
tained by  use  of  the  field  equations  and  the  series  is  telescoped;  3)  applica- 
tion of  interface  continuity  conditions  lead  to  differential  (in  time)  - 
finite  difference  (in  space)  equations;  4)  an  associated  set  of  differential 
tin  time)  - functional  difference  equations  is  observed  to  contain,  as  a 
solution  subset,  solutions  of  the  finite  difference  relations;  5)  solutions 

of  the  functional  equations  are  smooth  fields  which  take  on  exact  values 
at  discrete  spatial  points;  these  smooth  fields  are  now  employed  as 
dependent  variables;  6)  all  spatial  differences  are  expanded  in  a Taylor 
series  - an  operation  which  furnishes  partial  differential  equations,  i.  e.  , 
a continuum  model;  7)  nondimensional  variables  and  a small  parameter 
representing  the  ratio  of  typical  micro-to-macro  dimensions  is  introduced 
and  the  operators  are  truncated  according  to  an  asymptotic  procedure; 

S)  mixture  equations  result  by  algebraic  manipulation. 

The  second  procedure,  illustrated  in  Section  4 by  an  elementary- 
example,  is  the  Regular  Asyniptotic  Method.  This  approach  has  been 
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used  to  model  nonlLnea  r wave  propagation  in  laminated  [ 16]  and  fibrous 
]17-19]  composites,  and  the  rmal  diffusion  in  fibrous  composites  [20-Zl]. 
The  method  consists  of  six  basic  construction  steps  which  may  be  sum- 
marized as  follows:  1)  based  upon  the  particular  composite  and  the 
process  considered,  an  estimate  of  the  order  of  magnitude  of  all  field 
variables  is  made;  2)  the  field  equations  are  re-scaled  such  that  all 
dependent  variables  are  0(1);  in  the  process  a smiall  parameter  is  intro- 
duced which  constitutes  a measure  of  the  ratio  of  micro- to-macro- 
dimensions  of  the  problem;  3)  the  conservation  equations  are  averaged 
to  obtain  a standard  mixture  form;  4)  the  scaled  field  variables  are 
expanded  in  a regular  asymptotic  series;  5)  the  two  lowest  order  systems 
(including  interface  conditions)  are  used  to  obtain  expressions  for  the 
interaction  terms  and  the  constitutive  relations;  6)  the  constants  involved 
in  5)  are  obtained  from  the  solution  of  a static  micro  boundary  value 
problem  defined  over  a unit  cell;  a variational  principle  - based  finite 
element  method  is  suggested  to  resolve  this  boundary  value  problem. 


A judicious  choice  of  a particular  method  depends  upon  the 
problem.  The  Differential-difference  Method  is  best  suited  to  linear 
problcmj  and  elementary  geometries  such  as  lam.inates.  Here  formally 
exact  relations  may  be  obtained  in  many  cases.  The  Regular  Asymptotic 
Method  is,  on  the  other  hand,  applicable  to  nonlinear  problems  and 
complex  geometries.  Ln  this  case  one  is  usually  content  with  the  first 
few  terms  of  the  asymptotic  series. 

3.  EXAMPLE  I:  WAVE  PROPAGATION 

For  the  purpose  of  illustrating  the  Diffe rential- diffe rence  Method, 
consider  wave  propagation  in  a periodic  array  of  linearly  elastic,  homo- 
geneous and  isotropic  laminae,  perfectly  bonded  at  all  interfaces  as 
illustrated  in  Fig.  1.  Let  y = x^  for  notational  convenience  and  assume 
a state  o:  uniaxral  motion  in  the  y-direction;  i.  e.  , the  case  of  one- 
dimensional longitudinal  waves  propagating  normal  to  the  laminae.  Re- 
stricting attention  to  small  isothermal  deformations,  the  appropriate 
conservation  and  constitutive  equations  are 

= 0 , (a-{>.-2p)»^u)^“'^^  = 0 , (1) 

while  the  interface  continuity  conditions  are 


[u,cl  (h  ,t) 


th  th 

In  the  above  the  superscripts  ((i,k)  refer  to  the  k layer  of  the  Q 
constituent,  a = 1 , 2 ; k = 1 , 2,  3.  . . . The  variable  yh3t>i^)  ^ local 

coordinate  with  origin  at  the  midplane  of  the  (a,k)  layer;  0 = O^Z 
u = U2  are  normal  stress  and  displacement  (subscripts  on  C and  u have 
been  dropped  for  brevity).  The  quantities  p(0(),  h(a),  are  mass 

density,  layer  half  thickness  and  Lame'  constants  for  the  a - constituent. 
In  addition  >y{  ) = >(  )/^y,  S^(  ) s d(  )/>t. 


Wave  motion  in  the  composite  is  completely  specified  by  (1),  (2) 
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together  with  initial  conditions  at  t=0  and  boundary  data  on  y~0,  H(H  s ®). 
The  objective  of  the  subsequent  work  is  to  replace  this  set  of  equations 
(infinite  if  H = =)  by  a single  set  of  partial  differential  equations  which 
represent  a continuum  mixture  theory. 

3. 1 Expansion,  Recurrence  Relations 

Let  us  expand  the  stress  and  displacement  about  each  lamina  mid- 
plane as  follows: 

= I g In'.  . (3) 

n=0 

{ ^ ^ ^ ) 

Here  g represents  any  of  cT  ’ , u . Upon  substituting  (3)  into  (I)  and 

equating  powers  of  one  obtains  the  differential  - recurrence 

relations 

' • %)  Vl)  • 


(a)  ~(a),  (a)  ;:{a)  (a) 

c = E /p  , E = (X+2p)  , (4b) 

whe  re  n = 0,  1,2...;  k = 1,2,3,..,;  a = 1,2.  With  use  of  the  recurrence 
relations  (4a),  the  series  (3)  may  be  telescoped  as  follows 

{u,aj  (y  , t)  = C 

(a.k)  (a,k)  1 2 (a,k) 

Ty  S {E  pa.u}  (5) 

where  C , S ’ are  formal  differential  operators  defined  by 

jQ.k)  u(  ^ (0:.k)^  , 

C s cosh(c  y a^)  , 

<,la.k)  (q)  * (o.k)  -I  . (a)  ^ (a,k) 

S £ (c  y a^)  sinh(c  y • (6) 

The  tunctions  CT(o/®’^^(t),  represent  stress  and  displacement  at 

the  midplane  of  the  (o!,k)  layer;  these  serve  as  dependent  variables  in  the 
ensuing  analysis. 


3.2  Transition  to  a Continuum  Theory 

Substitution  of  the  expansions  (5)  into  the  interface  conditions  (2) 
furnishes  the  following  differential  (in  time)  - finite  difference  (in  space) 
equations  with  dependent  variables  defined  at  discrete  points  along  the 
y-axis: 

Jcy),  fci.k+l)  (c,,k)  ,JP)(P.k)  Jq)~(:^)’'  (a),  (a,k+l)  (o.k) 

Ja)  (a,k+I)  (a.k)  (P)  (p,k)  (al  (alja)  2 (a,k+ 1)  (a.k) 

'“(01  *“(0)  ' - “(0)  - P - (^0)  1 = 0. 
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Q,  P = 1, 2;  Of  P. 

[£  one  now  defines  the  fields 

(y.t)  H 1^,0}  (y»t)  for  ye 


ip.k)  Jp)  JP,k),JP), 


= {0.0}  for  yq(Y 
O,  p = 1,2:  C:  ? p;  and  if  one  notes  that 

[u.CjjQ^  (t)  = {u.a}  {Y  , t) 
then  (7)  can  be  rewritten  in  the  form 


C i_u  (y+i,t)  + u (y-i,t)]-2C  u (y.t) 

-1 


, Y"  '+h''^')  , 


Y(2.M._^I,k^ 


(8) 

(9) 


,(a)~(a)  (al  (a),  , , ..  (a),  . „ 

- h E S ta  (y+i,t)  - a (y-£^.t)j  = 0 


Ja).  (a),  . ^ (a),  . ,_{P)  (P).  .> 

C [o  (y+^.t)  o (y-A.t)]  - 2C  a (y.t) 


(a)  (a)  (a)  2 (a)  (a) 

-p  h S >^[u  (y+A,t)  - u (y-£,,t)]  =0 


(10a) 


(10b) 


a.P  1.2;  a / p.  Equations  ( 10a,  b)  were  obtained  by  replacing  the  fields 
(-3),y(P,k)^  t)  and  + A.  t),  etc.,  by  the  smooth  fields  o(P)(y,  t), 

'^^Hv'A.t),  etc.,  defined  over  the  entire  domain  of  the  composite;  here 
■>.-(a,.<)  denotes  the  midplane  of  the  (a,k)  layer  (see  Fig.  I). 

The  functio nal  difference  equations  (I0a,b)  contain,  as  a subset, 
the  solutions  of  the  finite  difference  equations  (7).  That  is,  the  new 
lields  satisfy  (7)  at  y = 'i'tP.k).  Therefore,  the  new  dependent  variables 
if^l,  c(a)  are  functions  defined  for  all  y which  assume  exact  values  at 
layer  midolanes.  Assuming  the  admissibility  of  such  expansions,  all 
s’urr.s  and  differences  are  now  expanded  about  A ==  0 in  a Taylor  series; 
Enis  leads  to  the  partial  differential  equations 

,(0)^(0)  Jo).  -(ol'^Jo)-  (a)  Jp)  (P)  ^ 

-Cu  =0, 

Y 


Cu 


h'  J E 


Jo)-  (a)  (o)  (o)  (o)-  2 (o) 

C Ca  - h ap  S S>  ? u 

Y t 


JP)  (P)  . 

c a =0 


(11) 


a,  P = 1,  2;  a ^ P,  where  C,  S are  formal  differential  operators  of  the  form 


C = cosh(A3  ) . S = (AS  ) ^ sinh(AS  ) • 

Y Y Y 

Equations  (II)  constitute  a continuum  model  of  the  composite.  The 
dependent  variables  u*^'*,  o^.o)  of  ( I 1 ) can  be  used  to  reconstruct  the 
microstructural  fields  through  (5)(noting  (8)). 

3.3  Scaling 


(12) 


Let  the  typical  ’ macrodimens  ion"  of  the  problem  be  A and  the 
typical  microdimension  be  A.  The  quantity  / represent  a dominant 
signal  wavelength  associated  with  wave  motion.  Introduce  the  non- 
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dimeasional  quantities 

Z = y/  ^ . T H tc/ / 


Ja)  (a).  2 (a) 

fsA/-'-  . I =0  /pc  . v 


(a). 
u 11. 


(13) 


where  p,  c denote  appropriate  mLxtore  density  and  wave  speed.  Then  (11) 
be  corne  s 


where 


Ja)  (a) 

C ==  cosh(ry  > ^) 


Ja)  ,,  (a).  >- 1 . , , (a)  . 

5 = (cy  S^)  sinh(fy  ;^^) 


C £ cosh(fa^)  , i = (ca^)  ^ sinh(fa^)  , = (c/c'^^n*^^ 

(O')  ,(a),  —(a)  (a)  (a),  -(a)  2 (o;)  ~{a) 

n = h /A  , P = n p /p  , J s pc  n /E 


(14) 

(15a) 

(15b) 

(15c) 


The  ?overnin5  equations  ('4)  are  now  scaled  such  that,  if  y^(0,/) 
and  Tc(0.  1/ c),  then  ^(^(0,  I)  and  Tc(0,  I);  that  is,  the  macrodimension  is 
now  0(1)  while  the  microdimension  is  0(f).  Thus,  it  is  evident  that  all 
roregoing  Taylor  expansions  about  *1=0  are  in  fact  expansions  about  € = 0; 
further,  all  operators  are  power  series  n =0,1,2,  ..  . Such 

expansions  are  not  required  to  converge,  but  only  be  asymptotic  as 
r -*  0 *'4'. 


3.4  Binary  Mixture  Theories 

With  some  algebraic  manipulations,  equations  (14)  may  be  re- 
written in  a form  characteristic  of  binary  - type  mixture  theories  as 
follows: 


, In)  Ja)  .^In)  (al-  2 la).  , , a+ 1 

C - p "S  . = (- 1)  S , 

C n -e  S SL  .=(-1)  p. 


= 1,2.  He  re 


(16a) 

(16b) 


C"’  =1/21  + (fB^)  /41  + (fB  )*^/61...  , 


(16c, d) 
(16e) 


The  bracket  in  (16a)  suggests  a momentum  equation  while  that  of 
(I  6b)  has  the  character  of  a constitutive  relation.  The  quantities  p and  R 
represent  constituent  interaction  terms. 


A m.ixture  theory  of  order  n is  now  defined  as  those  equations 
obtained  from  (16a,  b)  by  truncating  the  operators  ' and  as 

well  as  the  in  p and  R,  after  0(f“-’^‘’)  terms.  In  particular,  a 
rr.ixture  theory  of  order  1 is,  under  this  definition,  given  by 
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a (?  r/2  - = (-1)“^^'= , ^A^^nf2-pj:{'‘^ 

K C 'P  ^ ^ 


{-  1 ( 17a, b) 


2 

c“p- 


(18c. d) 


Equations  (17a)  may  be  interpreted  as  mixture  momentum  relations  and 
(17b)  as  mixture  constitutive  equations.  The  interaction  terms  in  this 
formulation  are  proportional  to  displacement  and  stress  differences, 
respectively. 


The  operator  preceding  the  left  side  of  each  equation  (17a, b) 
can  be  eliminated  by  expanding  the  dependent  variables  in  regular 
asymptotic  series  in  f.  The  result.upon  retaining  O(f^)  terms  only,  is 


,(“P) 


(ap)  2 (a) 

- P 


(18) 


where  jg  is  a constant  which  is  a complicated  function  of  composite 
geom.etry  and  material  properties.  The  quantities  p(<^P)  here  play 

the  role  of  "partial"  stresses  and  densities;  their  physical  interpretation, 
however,  is  not  obvious.  Details  concerning  the  derivation  of  (18), 
which  shall  be  called  a modified  (first  order)  mixture  theory,  can  he 
found  in  [41. 


Let  us  now  consider  the  limit  as  f -♦  0 in  (16).  If  it  is  assumed, 
that  t and  c may  be  selected  such  that  o^(  )=0(1),  S»(  )=0(1), 

n = 1,2...,  whe  re  ( ) represents  t/®  ' or  as  f-»0(i.e.,  as  the  ratio 

of  micro-to  macro-dimensions  vanishes),  then  (16)  reduces  in  dimen- 
s ional  fo  rm  to 


-2  , (1)  (1),  (2)  (2)^  (I)~(l) 

IS  - (n  p rn  p )(n  E 


- 1 


+ n 


(2)~(2)'\  2, 


,(o) 


)?,]{u.ar  = 0.  (19) 


Therefore,  the  theory  is  nondispe  r sive  as  c -»  0,  a result  to  be  expected. 
Equation  (19)  implies  that  the  appropriate  mixture  density  p,  modulus 
E,  and  wave  speed  c,  are 

(1)  (1)^  (2)  (2)  ~-l  (1)~(1)“  (2)~(2)'^  2~ 

p = n p +n  p , E = n t.  +n  E , c =E/p. 


(20) 


3.5  Phase  Velocity  Spectra 

It  is  evident  that  the  foregoing  theory  represents  a long  wave- 
length, low  frequency  approximation  of  composite  wave  motion.  A 
natural  question  concerns  the  range  of  validity  of  this  approximation.  A 
convenient  indication  of  the  comparative  accuracy  of  an  approximate 
theory  can  be  obtained  by  studying  the  phase  velocity  spectrum. 

Assuming  steady-state  sinusodal  wave  trains,  equations  (16),  (17)  or  (18) 
become  a set  of  homogeneous  algebraic  equations  from  which  the  relation 
between  wave  number  k and  phase  velocity  Cp  is  calculated.  A typical 
comparison  of  exact,  first  order,  and  modified  (first  order)  mixture  phase 
velocity  spectra  is  shown  in  Fig.  2 for  a Polymethyl  methacrylate 
(PMMA)  - 304  stainless  steel  composite.  The  data  shown  represents 
the  first  n\ode  of  propagation  down  to  the  cutoff  frequency  associated 
with  the  first  "stopband  ’.  The  exact  spectrum  results  from  retaining 
all  te  rms  in  the  ope  rato  rs  of  (16)  and  is  identical  to  that  g iven  by 
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Ry-tov  [22].  Comparisoas  such  as  Fig.  2 indicate  that  reasonable 
accuracy  may  be  expected  from  a theory  of  order  I down  to  wavelengths 
equal  to  several  microdimensions  for  problems  where  the  energy  is 
partitioned  primarily  in  the  first  mode.  This  is  quite  a remarkable 
result  for  such  an  elementary  theory.  Additional  details  concerning  the 
phase  velocity  spectrum  as  well  as  transient  pulse  calculations  and  an 
associated  discussion  concerning  initial  and  boundary  conditions- can  be 
round  in  ;_4-9]. 

4.  EXAMPLE  2:  HEAT  CONDUCTION 

^ For  the  purpose  of  illustrating  the  Regular  Asymptotic  Method, 
consider  a linear  thermal  diffusion  process  in  a periodic,  two-dimensional 
array  of  unidirectional  fibers  of  arbitrary  cross  section  embedded  in  a 
matrix,  as  Ulustrated  in  Fig.  3a.  Let  a ''cell”  be  associated  with  each 
tibe^r  as  depicted  in  Fig.  3b.  Each  such  cell  consists  of  regions  A<^^  and 
occupied  by  the  fiber  and  matrix,  respectively,  with  unit  outer 
normals  N^  > and  N^^).  The  interface  between  the  two  constituents 
shall  be  denoted  by  1,  and  the  outer  boundary  of  the  cell  by  C.  With 
respect  to  rectangular  Cartesian  coordinates  xi,  x^,  xj  as  shown  in 
r >.5.  j<i.  icL  me  coiiipusiie  occupy  the  domain  Ci  s.  x^  S L,  - =,  < < m , 

- = < xj  < = . Assume  that  the  two  constituents  in  each  cell  are  homo- 
geneous and  isotropic  with  a perfect  interface  (no  thermal  resistance). 
Finally,  let  the  initial  conditions,  and  the  boundary  conditions  on 
X3  - 0,  L be  such  tliat  the  temperature  field  is  similar  in  each  cell.  In 
view  of  the  lastpremise  it  is  sufficient  to  consider  a typical  cell  with 
rero  heat  tlux  normal  to  the  boundary  C.  Consequently  the  basic  equations 
lor  the  temperature  fields  T^“'  and  heat  flux  vectors  Q-(=^).  o = 1 2 are 


(Q.  . + dt)'^^  = 0 
J.  J 


(Q.  + kT,  = 0 

i i 


on  A 


(o) 


(21a.b) 


(q'^'n'^Lo 

P i3 


on  C:  f 


- T and  (Q. 

|3 


(2)  (!’ 
s H ' 


0 on  I. 


Equation  (21a)  is  the  conservations  of  energy;  the  relation  (21b)  is  thl  ’ ^ 
rourier  heat  conduction  law;  equation  (2lc)  is  a symmetry  condition;  and 
(-Id)  represents  tne  interface  continuity  of  temperature  and  normal  heat 
£*ux  (note  N3>  = G).  The  superscript  a = 1,2  refers  to  the  fiber  and  the 

m.atrix,  respectively.  The  quantities  d,  k denote  heat  capacity  and 
tnermal  conductivity,  respectively.  The  notations!  ),:=»{  )/ax  and 
( ) = 3(  )/ St  have  been  introduced  for  convenience. 

4. 1 Seal ino 


Equations  (2  1),  together  with  initial  conditions  at  t=0  and  appropri- 
ate boundary  data  on  xj  = 0,  L specify  a well  posed  problem  involving 
tnr^  spatial  variables  x,  and  time.  t.  The  objective  of  the  subsequent 
analysis  IS  to  derive  a continuum  mixture  theory  governing  the  macro- 
^copic  diffusion  process  which  involves  only  one  spatial  dimension  (X3) 


-Latin  and  Greek  subscripts  imply  Cartesian  tensors  and  the  usual  sum- 
mation convention  witn  ranges  from  1 to  3 and  1 to  2 respectively. 
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and  vet  reflects,  at  least  approximately,  the  effect  of  conduction  on  the 
microscale.  To  this  end,  let  A and  A be  associated  with  typical 
mac  roscopic”  and  'microscopic''  observational  dimensions,  respectively-. 
These  lengths  m.ay  be  definea  in  terms  of  characteristic  thermal  diffusion 
tim.es  in  the  longitudinal  and  transverse  directions  according  to 


t,  . , H d .A  /k 
(A)  (m)  (m) 


2 

(L)  (m)  (m) 


— (22) 


where  d(ni}<  ^(m)  mixture  heat  capacity  and  thermal  conductivity 

(to  be  defined  later).  In  addition,  it  will  be  convenient  to  introduce  a 
mixture  heat  flux  Q(nT.)  based  upon  a reference  temperature  T,  and  a 
narameter  f as  follows: 


Q,  , s k T/A 

(m)  (m) 


(23) 


The  quantity  f represents  the  ratio  of  micro-to-macro-dimensions  of  the 
com.pos  ite. 


With  the  aid  of  the  foregoing  definitions,  the  following  non- 
dirre  ns  ional  variables  a.re  introduced 


T sT  /T,  : sd  /d,  ,,  ft  sk  /k, 

m m 


(23a) 

(2.3b) 


whence  (21)  become 


• (o ) 2 (o()  (a)  (a)  (o) 

(D.  +iT)'  =0  , (-  . c n )'  ' + ft'“'(7,  = 0 ontf  . (2.4a, b) 

J»  J 3 P 3 p 


n<“'  o n J2)  (1)  (2)  (1) 

2,  V = 0 on  C : o = o . (T,  -Do  )^r,  = 0 on  J . 

P P p p P 


(24c, P) 


Jy)  (o) 

:24)  P = I,  Z^and  C , J,  C denote  .-\  , 1,  C in  nondimens ional  form; 

the  vectors  are  unit  outer  normals  to  the  boundaries  of  Ci^^;  partial 

deri-.-atives  are  now  defined  by  ( ),  , s §(  )/d^.,  (*)  s b(  )/2T  . 

4. 2 MLxture  Equations 

Binary  mLxture  equations  for  this  problem  can  be  immediately 
obtained  by  averaging  (24a)  overG^'^f)  according  to 

- I 

(C.)  (o)  (o) 

wnere  3 represents  D or  G • This  yields,  after  a little  algebraic 
m.anipulation 


3.  3 


)3  (3 


= G ^ D , V ds  , 

P P 


(26a,b) 
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whe 


(ctp)  (a)^(aa) 
-3  = ^ D3 


(ap)  (a)  (a) 
d = n ^ 


n £ G /(G  +C  ) . 


(ap ) 

Tae  quantities  2 3 


(ap)  (a) 


n represent  partial  heat  fluxes,  partial 
heat  capacities,  and  volume  fractions,  respectively.  The  variable  p is 
an  interaction  term  reflecting  heat  transfer  from  the  fiber  to  the  matrix. 
At  this  stage  the  mixture  equations  (26)  are  exact.  Approximations  arise 
when  one  attempts  to  m.odel  the  Fourier  expressions  for  heat  conduction, 
and  the  interaction  term. 

4. 3 Asymptotic  Expansions 

A fundamental  premise  is  now  introduced:  the  ratio  of  the 
characteristic  thermal  diffusion  times  in  the  transverse  and  longitudinal 
(fiber  axis)  directions  in  small  compared  to  unity,  i.  e.  , 


(28) 


Equation  (28)  is  appropriate  for  many  composites  used  for  thermal  pro- 
tection. The  premise  (23)  suggests  the  following  regular  asymptotic 
expansion  for  all  dependent  variables,  denoted  by  G'*-''': 


(a) 

d '(C..T;f) 


^ 2n  (a)  , 

L f G/1^,t)  . 

n 2n  ^1 

n=0 


(29) 


If  (29)  is  substituted  into  the  governing  equations  (24a,  b)  and  the  coef- 
ficients of  similar  f-order  are  equated,  one  obtains  a system  of  equations 
for  each  n = 0,  I,  2.  . . . In  what  follows,  a mLxture  theory  is  developed 

based  on  the  lowest  order  system. 

4.4  Interaction  Term 

The  lowest  order  system  corresponding  to  (24a,  b)  is 


(2. 


12- 


3\0) 


(c^ ) 

(T.  ^ - G, 2 - °\o) 


‘J.  J ' ' (0) 

Equations  (30c)  yield 

„(a)  _ ^(a), 

'^(0)  ^ ‘^(0)'^3'  ’ 

whence,  to  lowest  order  accuracy,  equations  (25),  (27)  furnish 


(30a, b,c) 


(31) 


= - * 


(ap)^  (ap) 

3 '-^3 


O . (^  -,  » T 


(32) 


In  order  to  complete  the  mixture  formulation,  the  functional 
dependence  of  the  interaction  term  P on  the  averaged  temperatures 


Jaa) 

;'(a) 


must  be  determined.  For  this  purpose  it  is  necessary  to  consider 

and  to  satisfy  the  continuity  of  temperature  including  O(f^)  terms. 
To  begin  this  task,  one  finds  from  (32)  and  (24b)  that 


^fa) 

13(0) 


^(a)  (a) 

• ~{2),3 


(33) 


Next,  with  use  of  (33)  and  (30b,  c),  equation  (30a)  becomes 
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(a)  (a) 

^ • <34) 

(q) 

i ne  functions  o can  be  related  to  n,  defined  by  (26b),  by  integration  of 
'.34)  ov'er  G*  K use  of  Gauss'  Theorem.,  and  the  syinmetry  condition  (24c); 
the  result  is: 

(a)  (a)  , , ,a 

n o = (-  1)  P . (35) 

Thus,  one  obtains  the  following  relations  for 

^ (2) 

(a)  (a)  (a)  . , , ,a  (a)  „ 

° * ^(2)’p  'p  ' ^ P on  G (a,P  =1.2)  . (36) 

1 he  appropriate  boundary  conditions  are,  from  (24c.  d). 


3'^’  ,v'2>=0  on  C,  (."U"  = 

(2)  p p ^ ■ ■"(2). p • '-'(2). a a 


~(2)  P P - "(2).p  '■'(2).p'‘p  ’ 

J2)  (1)  _ .Jl)  _(2)^,  2 

^(2)  ^(2)  '"*(0)  ■ "(0)  on  J . (37) 

£,quaciuns  loo^,  \yii  constitute  a nucto'  boundary  value  problem 
^7^  ^1’  ^2  ” plane;  the  solution  is  unique  within  a function  V(^2,  T).Once 

are  known,  one  can  write 


A,-.he  re 


(al  (a)  2 ~(a)  2 (a)  4 

7 = ■ s-Q)  + f Jf  = ~(0)(^3’  ^ f r.7-  '(4]’42^  ^ ^ 

p-4=<‘>=J<>  n--<^>  - (2)  ^ ,-(2)  ^(1)  2 

■ ‘^(2)  ’ = ^(2)  ■ '^(0)  ■ "(0)^'  f • 


Upon  averaging  (39),  one  obtains 

1 , 2 ,.(aa) 

a r)  - t)  + ( P7-  . (40) 

r;ence,  the  interaction  term  can  be  written 

p _ ,,JIa)  J2a)  2 ,„Jla)  J2a)-I 

P - a{a  -o  )/f  , a = (G-  -G*  ) . (41a,  b) 

Equations  (41a),  (32),  (26a)  close  the  mixture  theory.  The  primary 
result  is  expression  (41b),  which  can  be  used  to  determine  the  interaction 
coefficient,  a,  for  arbitrary  fiber  and/or  cell  geometry  via  the  solution 
of  a time  independent  problem  defined  over  the  unit  cell.  In  terms  of 
7'*  ^ . this  problem  takes  the  form; 

3f' o„  r : 3 = = 0 o„.,  H„ 

7"'^ ' *=0  at  point  Ocr  * * . (44) 

In  general  it  is  difficult  to  obtain  an  analytical  solution  to  the 
above  micro  boundary  value  problem.  As  an  alternative,  a finite 
element  procedure  has  been  proposed  /2I]  which  may  be  utilized  for 
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arbitrary  two-dimensional  cell  geometry.  This  procedure  Ls  based  upon 
a modified  Re  iss  ne  r -type  variational  principle,  the  corresponding 
functional  of  which  is  defined  by 


a = l (of) 

n 


o-o  + — 7—. — jdG  + oD''(~-''  -T-'  )ds 

a J (45) 


{ ^ } ^ cy ) 

vshich  is  defined  in  terms  of  on  G and  2^-  on  J.  The  equations  for 

-.iof)  £oiiow  from  the  requirement  H be  stationary  with  respect  to 
arbitrary  variations  of  and  2";  here  2"'  is  a Lagrange  multiplier 

which  physically  represents  the  heat  flux  normal  to  the  boundary  J. 


4.5  Numerical  Results 

The  finite  element  procedure  has  been  used  for  the  solution  of  the 
temperature  microstructure  problem  and  the  interaction  coefficient  for  a 
variety  of  geometries  and  combination  of  material  properties.  Details  of 
such  calculations  can  be  found  in  [21].  Here  several  results  of  particular 
interest  are  noted;  properties  used  are  listed  in  Table  3. 

A frequent  fibrous  composite  geometry  is  a hexagonal  array  of 
circular  cylindrical  fibers.  In  practice  it  is  common  to  approximate  the 
cell  geom.etry  of  this  case  by  concentric  circular  cylinders  [20]. 
N’umerical  results,  based  on  the  theory  presented  herein,  related  to  the 
accuracy  of  this  approximation  are  shown  in  Fig.  4,  where  the  inter- 
action coefficient  has  been  calculated  as  a function  of  fiber  volume 
fraction  and  thermal  conductivity.  From  the  data,  the  conclusion  can  be 
drawn  that  for  a practical  range  of  fiber  volume  fra ctions , the  concentric 
cylinder  approximation  can  be  used  without  any  significant  loss  of  ac- 
curacy. 

For  calculation  of  the  interaction  coefficient,  the  concentric 
circular  cylinders  approximation  based  on  equal  fiber  volume  fraction 
can  also  be  used  for  composites  containing  square  fibers  arranged  in  a 
square  array.  This  is  borne  out  by  the  results  shown  in  Fig.  5. 

Aa} 

The  fields  o'*'  (4l’C2^  needed  to  calculate  the  interaction  term  can 

also  be  used  for  the  calculation  of  the  temperature  distribution  in  a unit 
cell.  From  (38)-(40)  one  has 


7 (^■. 


. Jna) 
t)  - o 


Ln  the  right  side  of  (46)  the  only  quantity  dependent  upon  inplane  coor- 
dinates ^1  and  ^2  It  follows,  therefore,  that  curves  of  equal 

2 are  also  isotherm.al  lines  within  the  fram.ework  of  the  mixture 
theory.  For  this  reason,  and  to  illustrate  the  type  of  temperature  micro- 
structure  that  can  be  obtained  from  the  mixture  theory,  contours  of  equal 
r suitably  normalized,  are  given  in  Fig.  6 for  square  fibers  in  a 
f imila r array. 


Finally,  an  indication  of  the  accuracy  of  the  proposed  mixture 
model  and  the  construction  procedure  can  be  obtained  by  comparison  of 
mixture-predicted  response  with  2D  or  3D  computer  code  calculations. 
This  has  been  accomplished  in  [20]  for  the  case  of  concentric  circular 
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cylinders  where  two  spatial  diifiensLons  and  time  are  involved. 
Calculations  were  conducted  to  determine  the  evolution  of  the  temperature 
field  in  a quescent  half-space  s 0 subject  to  the  boundary  condition 


T{0,  r,t) 


1,  0 s t 5 t 

o 

0,  t > t 


(47) 


' o ' 

Both  the  mixture  equations  and  tlie  full  2D  problem  were  solved  by  finite 
difference  schemes  [20];  the  latter  was  used  as  the  correlating  norm  for 
estimating  the  accuracy  of  the  mixture  theory  solutions.  Fig.  7 illustrates 
typical  average  temperature  profiles  in  the  two  constituents  after  a short 
time  following  the  termination  of  the  temperature  pulse.  The  ability  of 
the  mixture  theory  to  predict  the  temperature  microstructure  is 
illustrated  in  Fig.  8;  the  radial  distribution  of  temperature  obtained  from 
the  mixture  theory  is  almost  identical  to  the  exact  solution. 


4.6  Remarks 

The  s i :"'nifica  nt  feature  of  the  above  mixture  theory,  and  of  similar 
theories  proposed  for  wave  propagation  [16-19],  is  that  it  converts  what 
is  essentially  a three-dimensional  problem  to  a problem  involving  a 
single  spatial  variable,  without  losing  the  essential  details  of  local  field 
distributions.  The  reduction  in  the  number  of  dependent  variables  leads, 
of  course,  to  a substantial  increase  in  numerical  efficiency. 


5.  SELECTED  APPLICATIONS 

Two  applications  involving  nonlinear  wave  propagation  are  briefly 
reviewed  in  this  section  in  an  effort  to  demonstrate  the  utility  of  mixture 
theories  with  microstructure.  The  first  involv-es  delamination  of 
iam.inated  composites.  The  second  concerns  finite  amplitude  elastic- 
plastic  wave  propagation  in  fibrous  composites. 


5. 1 Wave  Propagation  and  Debonding  in  Lam.inated  Composites 

A "first  order"  mixture  theory  with  microstructure  has  been 
developed  in  [12]  for  longitudinal  wave  propagation  and  debonding  in  a 
laminated,  binary  composite  with  periodic  microstructure.  The  case 
considered  concerns  small  isothermal  deformations,  linear  elastic  con- 
stituents, and  a Mohr-Coulomb  interface  failure  and  slip  model.  In  view 
of  the  last  item  the  problem  is  nonlinear.  With  respect  to  Fig.  1,  a 
condition  of  plane  strain  was  assumed  in  the  X3-direction  and  wave 
motion  yielding  symmetric  velocity  v^  and  anti-symmetric  velocity 
distributions  with  respect  to  the  X2  - coo  r dinate  within  each  lamina. 
Theoretical  construction  was  based  on  the  Differential-difference  Method. 
Ln  dimensional  form,  and  with  reference  to  Fig.  1,  the  relevant  equations 
a re 


L 
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(ap)  (orp)-  ,.<^+1 

M-I)  P 


’I 


(orp)  2 (q,  ) ( 

cr,,  = Ec  ^u/ 
y=l 


‘l.I 


(48a, b) 


P = K(u 


(la)  .(2al 


P = (■ 


1 1 
A-Ba.^ 


if  P<lP 


cr  ' 


■33gn(u^^^^-u|^^*)  if  P = 1 P 

11  ' cr 


(48c) 


Equations  (4Sa)  are  mixture  momentum  equations  while  (48b)  represent 

mixture  constitutive  relations.  The  quantity  P denotes  a momentum 

interaction  term.  The  variables  ai/“P>,  p(ap)  denote  partial  quantities  and 

are  defined  by  ( s where  ( )*“^^  is  an  averaged  variable 

and  a=^  = h ^Mh  > + The  quantity  1 P^j.  ] is  related  to  the  interface 

shear  stress  CTf->  and  normal  stress  O’’'  t>V 
12  22  ^ 


c r ' 


= cr 


I 7 


= A-Bo^2  ^''22  ^ ' 


rr-  - r (la)^  (2)  (2)^(1)  (2a),,  fl)  (2)  (2)  (I) 

*^22  * ' ^ ^ ^ ^ 1^^*"  ^ ^ > • <49) 

The  constants  A,  B represent  interface  cohesion  and  frictional  coefficients, 
respectively.  The  constants  K are  given  in  terms  of  the  composite 

properties  and  geometry  by 


i\pOu  ,1  /|h  -h  )(n  |i  +n  p ) 


(aar)  (Q')Ja)  X 
c s n E - — 


(o)' 


1.^. , _ E' 


,(1)  J2) 


E ’ “ 11) 

n' 


(2) 


, E^  ^ = (>+2p)^  \ (a,  P = I,2;  nf'P), 

(50) 


.■here  X,  P are  Lame  ' constants. 


Considerable  insight  into  the  modeling  capability  of  the  above 
mixture  theory  can  be  gained  by  comparison  of  theoretical  calculations 
with  the  experimental  Lm.pact  data  on  delaminated  plates  reported  in  [23]. 
The  composites  in  this  study  consisted  of  alternating  layers  of  Poly- 
methyl methacrylate  (PMMA.  Rohm  and  Haas  type  A)  0.762  mm  thick, 
and  6061-T6  alum.in’um.  0.792  mm  thick.  The  laminae  of  the  composite 
were  oriented  perpendicular  to  the  impact  plane,  and  struck  by  a 
projectile  from  a 10cm  bore  light  gas  gun.  An  aluminum  buffer  plate 
1.0cm  thick  was  placed  at  the  rear  of  the  composite  to  improve 
planarity  of  the  transm.itted  wave  front.  A transparent  Dynasil  100 
w-indow  followed  the  buffer,  and  a thin  mirror  was  vapor  deposited  at  the 
buffer  windo'w  interface.  The  motion  of  the  mirror  was  monitored  to 
within  rl.5  ^ lO'SrrjTi  i^y-  rneans  of  a displacement  interferometer.  The 
experimental  configuration  reproduced  from  [23]  is  depicted  in  Pig.  9 . 
The  shot  m.atrix  is  summarized  in  Table  1. 

The  mixture  theory,  consisting  of  (48)  was  coded  in  finite  dif- 
ference form,  and  solved  numerically.  Material  properties  for  the 
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composite  laminates,  flyer  plate  material,  buffer  and  window  materials 
are  listed  in  Table  I.  Input  to  the  computer  code  required  only  the 
material  properties  of  all  constituents,  goemetries  of  the  test  configura- 
tion, and  flyer  impact  velocity.  Numerical  values  for  A and  B were 
determined  by  a parametric  study  on  the  results  of  experimental  shot 
number  2.  They  were  found  to  be  A = O.OI  x lO'^  dynes/cm^  and  B = 0.50. 
All  other  correlations  were  performed  using  these  values.  Typical 
results  are  depicted  in  Figs.  10-12.  As  can  be  seen,  agreement  is 
generally  excellent.  Both  theoretical  and  experimental  results  clearly 
exhibit  the  effects  of  bond  breakage  and  delamination.  This  is  manifested 
by  the  appearance  of  a precursor  at  the  front  of  the  wave  profile  and  by 
greatly  reduced  oscillations  behind  the  wave  front  when  compared  with 
both  theoretical  and  experimental  results  for  perfectly  bonded  wave- 
guides (e.  g.  , see  [12]). 


5.2  Nonlinear  Wave  Propagation  in  Fibrous  Composites 

A binary,  first-order  mixture  theory  with  microstructure  has 
been  developed  in  [ 17],  using  the  Regular  Asymptotic  Method  and  the 
concentric  cylinders  approximation,  for  nonlinear  waveguide-type 
propagation  parallel  to  me  libers  ot  a unidirectionaliy  reinforced  fibrous 
composite  with  initial  periodic  hexagonal  array.  Fig.  13.  The  resulting 
m.odel  incorporates  the  effects  of  thermodynamics,  finite  deformation, 
and  elastic-plastic  constituents.  The  constitutive  relations  for  this 
analysis  consist  of  a Mie  - Gr  une  isen  equation  of  state  v/hich  relates  mean 
pressure,  density,  and  Internal  energy,  and  a von  Mises-type  yield 
--cirteria  and  associated  (J^)  floiv  rule  which  govern  the  stress  derivators. 
Lnder  the  assumption  of  an  adiabatic  process  and  isotropic,  homogeneous 
constituents,  the  following  mixture  theory  w'as  obtained  for  composites 
w-lth  perfect  bonds 

(a)  Continuity: 


la)  (aa)  iaa)  _ (aa)  (aa)  . • (aa)  , (aa) 

h ” i -U  , 1 = - (p  V'  ),  ^ T p ,/  p 


(51) 


(b)  Momentum: 


r (<^1  „.c  _(a)  (aa)  (aa)  (aa)_ 

[n  '-P-^S^^)  ].  3-n  P v^  =(-l)  P , 

whe  re 


^ (a)  Ja)  (aa) 

In  B £ 


r=l 


(a) 


(51a) 


(la)  (2a)3 

, ,,a  (a)Ja)Jaa)  2 ''^2  ^ 

- ( - 1)  n B f,  ^ = 


(51b) 

(51c) 
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ne^^emie  r 


Jl)  , J2)  ,,,  (2)^  3 (1)  1 (1)^  (1) 

B E I , R E - (2/n  )(—  - 2n  +—  n +ln  n ) 

(c)  Eaersv: 


2 (a(L  (aa)  (aa)  (aa)  (aa)  (ora)  , , >a+ 1 

t n [p  e - (-p+S^^  3]=(-U  R . 

0 = 1 

„ , c ,(la)  (1)  (la) 

R = (-pxS  ) n f 
rr 

(d)  Caloric  equation  of  state  (Mie- Gr^lnelsen); 

(aa)  (Q'a),,_,  (aa)  (aa)  (aa) 

P = g(p  )-rr(p  )p  e 

(e)  Stress  derivator  constitutive  relations: 

(aa)  (aa)  4 (aa)  1/aa)  Jla)  (la)  l_(la) 

33  3 3, 3 2 SS  rr  233 


(52a) 

(52b) 

(53) 

(54a) 


(2a)„{2a)  2 (aa)  1 (aa)  .3  1 

^ =--fv.  (W— In  — )1  . 


(54b) 


•-vhere 


9S  rr  33 


(aa)  ^(aa)  I (a)  (aa)  (aa) 

A = 0 if  ^ "3  ^ or  < 0 , 


(54c) 


Ja)  , (aa)  (aa)  (aa)  (aa)  _(aa)  (aa) 

^ ^33  "3,3"^r  ^ ^2 


and  J 


(aa  I Y 


.(a) 


2 3 

(f)  Interface  stress  nnatch: 


(54dj 


(-p-rS  ,t)  = (-p+S  ,t) 

rr  3 rr  3 

(g)  Volume  fractions: 


(55) 


^(la)  (1)  (l)dla)  (1)  (2)  , 

n = n f , n +n  = 1 . 


(56) 


(h)  Definition  of  averages: 


(1) 


( ) 


(la)  I 


TT  r 


(1) 


2 0 


2TTr(  )^  ^^dr  , 


(57a) 
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I ) 


,(2a) 


2.r(  /"’dr 


Ti(r 


(2) 


(1)  r (I) 

■r  ) r 


(57b) 


The  dependent  variables  above  number  14.  These  are  , 

and  e(°'^)(a  = l,  Z).  The  foregoing  represent  volume 
fraction, density,  velocity  component  in  the  X3-dir ection,  stress  deviator 
component  in  the  X3 -dire ction,  and  internal  energy,  respectively. 
Equations  (5I)-(56)  constitute  14  independent  relations  for  the  above 
unknov.'ns.  The  independent  variables,  which  were  originally  r,  B.x^,  t 
(see  Fig.  13),  are  now  x-j,t  only.  All  equations  are  written  in  an 
Eulerian  description.  It  should  be  noted  that  or  = 1 refers  to  the  fiber 
while  a=2  denotes  the  matrix. 


In  an  effort  to  test  the  validity  of  the  first-order  mixture  theory,a 
comparison  was  made  between  calculations  performed  using  the  quasi- 
one-dimensional  mixture  equations,  and  a formal  two-dimensional 
solution  ot  tne  concentric  cylinders  problem.  i he  mixture  solutions  were 
obtained  by  writing  (50-56)  in  finite  difference  form  and  solving 
simultaneously  for  the  14  dependent  variables.  The  necessary  numerical 
analysis  and  the  resulting  code  is  described  in  [24].  The  two- 
dimensional  solutions  were  generated  using  a well-known  Lagranglan 
finite  difference  computer  code  called  CRAM  [25].  A typical  cell  of  a 
quartz-phenolio  composite  occupying  the  half-space  X3  < 0 was  selected 
for  study.  For  computation  in  the  range  0 s p s 30  K bar  a Mie- 
Griineisen  relation  of  the  form 


(aa)  3 (a)  (2)  (a)  (aa)  (oa) 

P - IH  T?+Gp  e 

i=l  ^ 


(a) 


(ota),  (a)  , 

.p 


(58) 


v/as  used,  where  denotes  the  initial  density  of  mate  rial  o,  and 
G(^^  are  material  constants.  The  constituent  material  properties 
employed  in  all  calculations  are  given  in  Table  2.  Quiescent  initial 
conditions  were  assumed.  Boundary  conditions  consisted  of  a particle 
velocity  of  0,5cm/p  sec  of  1.0-p  sec  duration  applied  uniformly  to  both 
constituents  at  x^  = 0. 


A comparison  of  mixture  theory  and  two-dimensional  finite  dif- 
ference code  predictions  of  averaged  constituent  particle  velocity  as  a 
function  of  time  is  illustrated  in  Figs.  14,  15.  The  two-dimensional  code 
prediction,  which  resolves  the  particle  velocity  distributions  in  the 
radial  direction,  where  averaged  according  to  (57)  in  order  to  provide 
direct  comparisons  with  the  mixture  theory.  Results  are  shown  for  a 
propagation  ( X3-dire ction)  distance  of  0.5cm.  Agreement  between  the 
finite  difference  and  mixture  theory  calculations  is  judged  to  be  excellent. 
Finally,  it  is  noted  that  on  the  same  computer  (aUNIVAC  1108),  the  CRAM 
calculation  required  roughly  500  times  the  computer  time  needed  for  the 
equivalent  mixture  calculation. 
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Table  1. 


The  Experimental  Configurations. 


Com.posite 

thickness 

(cm) 

Projectile 

Experiment 

No. 

Material 

V elocity 
(cm/  ps  ec) 

TTiicknes  s 
(cm) 

o 

00 

to 

o 

,A.luminum 

0, 001355 

1.  634 

2 

0.  776 

Alum.inum 

0. 00  1420 

1.  631 

3 

0.  806 

Aluminum 

0,00 1 144 

1.  571 

4 

0.  S03 

PMMA 

0.0021  18 

0.  691 

5 

0.  305 

PMMA 

0. 003078 

0.  694 

6 

0.  812 

Tungsten 

Ca  rhide 

0. 00  1289 

0.  975 

7 

0.  810 

Aluminum 

0, 00  1330 

0.656 

8 

0.  809 

Aluminum 

0. 00  1066 

0.  164 

In  all  experiments  the  buffer  plate  was  1 cm  thick.  The  bond 
thicknesses  v.  ere  all  less  than  0.0002  cm  thick. 

* Experim.cnt  1 exhibited  excessive  flyer  tilt  and  Drumheller  and 
liUnder^an  considered  the  data  useless. 
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Table  2 


Material  Properties  used  in  TINC  and  CRAAI  Calculations 

Material  1 Material  2 


( (X  ) 

Initial  Density  pj  (g/cc) 

1.  85 

1. 29 

Shear  Modulus,  p.^^^(dynes/cm^) 

106. 0 (10®) 

30.  3 (10®) 

r ) 

Bulk  Xlodulus,  (dynes/cm^) 

397.  6(10®) 

96. 0 (10®) 

M J ^ ( dyn  e s / cm  ® ) 

38.0  (10^^) 

0 

^ e s / cm  ® ) 

27  5.0  (lO'M 

11.  3(10“) 

Mie-Griineisen  Ratio, 

. 32 

. 52 

Yield  Stress  in  Simple  Tension 

(dynes/cm^) 

2.75(10®) 

I.  33  (10®) 

Radial  r'iber  R.adius,  r^  (crti)  O.OS 

Outer  Cylinder  Radius  r (cml  0.  20 
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TABLE  3 


\IATERIAL  PROPERTIES  USED  FOR  COMPUTATIONS 


Thermal  Conduct! Ratio 


Specific  Heat  Ratio 


r.(l) 


50 


0.5 


\'olume  Fractions 


= .2  (Fiber),  = .8  (Matrix) 
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